The inner nebula and central binary of the symbiotic star HM Sge 
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ABSTRACT 

We present contemporaneous HST WFPC2 and VLA observations of the symbiotic nova 
HM Sge. We identify a number of discreet features at spatial scales smaller than ~ 0.1 arcsec 
embedded in the extended nebula, with radio and optical emission well correlated in the inner 
1 arcsec. For the first time we measure the positions of the binary components of a symbiotic star 
directly. We estimate the projected angular binary separation to be 40±9 milli-arcsec, with the 
binary axis at position angle 130°±10°. The latter is consistent with previous estimates made 
by indirect methods. The binary separation is consistent with a previous estimate of 50 au if the 
distance is 1250±280 pc. Temperature and density diagnostics show two distinct regions in the 
surrounding nebula, with a cool wedge to the south-west. An extinction map indicates the true 
interstellar extinction to be no more than E(B-V) = 0.35. This is consistent with a minimum 
distance of ~700 pc, but this would be reduced if there is a circumstellar contribution to the 
minimum in the extinction map. The extinction map also suggests a patchy dust distribution. 
We suggest that a southern concentration of dust and the south-west wedge are associated with 
the cool component wind. Alternatively, the southern dust concentration is the cause of the cool 
wedge, as it shields part of the nebula from the hot component radiation field. 

Subject headings: binaries: symbiotic - stars: individual (HM Sge) - circumstellar matter - radio con- 
tinuum: stars 

1. Introduction with an ionized component of the CC wind. An in- 

frared sub-classification has been made, dividing 
Symbiotic stars are an extreme case of interact- the clasg into D (usty)-types, and S(tellar)-types, 

ing binaries, with separations of a few to a few 10s based on the dommant contribution to the IR. 

of au. They comprise a cool component (CC), typ- Early mode i ing f the radio emission has shown 

ically a red giant or Mira-type variable, and a hot that D _ types typically have separations of 10 or 

component (HC) , usually a white dwarf associated more times that of S -types. Some of these objects 
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show very slow optical outbursts, similar to those 
of novae, but lasting decades (see e.g. Kenyon 
1986 and Mikolajewska 1997 for further discus- 
sion.) 

HM Sge is a D-type symbiotic star which un- 
derwent such an optical outburst in 1975. Since 
then, there has been evidence of dust obscuration 
events, as in 1985 (Munari & Whitelock 1989) 
when Ay ~ 13.5 (Whitelock 1988; Munari & 
Whitelock 1989), and in 1979-1980 when Ay ~ 12 
(Thronson & Harvey 1981). The different values 
of Ay have led to a suggestion of clumping in the 
dust (e.g. Thronson & Harvey 1981), which was 
also suggested by Richards et al. (1999) to explain 
the nebula structural changes seen in the radio. 
Various distance estimates to HM Sge are sum- 
marized by Ivison et al. (1991). Those authors' 
own estimate of 1300 pc was based on a reddening 
estimate of E(B-V) = 0.53. Kenyon, Fernandez- 
Castro & Stcnccl (1988) used far-infrared IRAS 
colors to derive the extinction and a distance of 
1800 pc. These estimates clearly depend on the 
effects of variable circumstellar extinction. Re- 
cently the Mira-type cool component has been 
clearly seen in spectro-polarimetry (Schmid et al. 
2000), indicating that dust obscuration has sub- 
sided somewhat. The spectra of Schmid et al. 
(2000) are dominated by the line emission, with a 
rising continuum beyond ^8000 A. They also sug- 
gested that the binary axis had a position angle of 
123° in 1998. This contradicts tentative conclu- 
sions by Richards et al. (1999) placing the Mira- 
type due north of the hot component (i.e. posi- 
tion angle ~ 0°). Many of the observations of the 
outburst and subsequent development of HM Sge 
are described in Nussbaumer & Vogel (1990) and 
Miirset & Nussbaumer (1994). 

Corradi et al. (1999) have conducted ground- 
based imaging of the extended nebula of HM Sge. 
They find a number of filaments and discrete fea- 
tures embedded in a nebula extending to ~ 13 arc- 
sec. This emission is extremely faint, and if the 
nebula is due to activity in the central binary, sug- 
gests a much longer history to such activity than 
observed since 1975. 

HM Sge was detected in the radio shortly af- 
ter the optical outburst (Fcldman 1977). Radio 
surveys showed HM Sge to be one of the bright- 
est radio emitters in the class (Seaquist, Taylor & 
Button 1984; Seaquist, Krogulec & Taylor 1993). 



Kwok, Bigncll & Purton (1984) found a diffuse 
halo ~0.5 arcsec in diameter, surrounding a cen- 
tral nebula ~0.15 arcsec in size. More recently, 
Eyres et al. (1995) demonstrated that the struc- 
ture seen in the radio could be correlated with 
HST images at the same spatial resolution, sug- 
gesting that at least the HC was centrally lo- 
cated in the inner nebula. Richards et al. (1999) 
traced the development of the radio structure over 
5 years, and related it to the binary motion. They 
derived a period of 80 years and a binary separa- 
tion of 50 au 1 . Non-thermal emission associated 
with an E-W outflow has also been found (Eyres 
et al. 1995), and Richards et al. (1999) showed that 
this is rare in symbiotic stars due to the relatively 
short timescale on which the emission dissipates. 

Here, we discuss our observations of the neb- 
ula associated with this star using WFPC2 on the 
Hubble Space Telescope (HST). We also present 
contemporaneous radio images of the inner neb- 
ula made with the Very Large Array (VLA). These 
observations allow us to determine the physical pa- 
rameters in the nebula, and to attempt to relate 
the conditions with the binary interaction and the 
outburst history of the system. 

2. Observations 

The HST observations were made on 1999 Oc- 
tober 22, as part of a GO program 8330 on sym- 
biotic stars, as shown in Table 1. Three orbits 
were allocated to HM Sge, and observations were 
made in seven filters, including F218W, F437N, 
F469N, F487N, F502N, F547M and F656N. The 
exposure times and dominant lines for images pre- 
sented here are given in Table 1. Further details 
of these filters are available from Biretta et al. 
(1996). The calibrated data were retrieved from 
the HST archive. Each exposure was executed in 
two sub-exposures to allow cosmic ray subtrac- 
tion. The F656N image presented here was also 
dithered to allow recovery of the full spatial reso- 
lution (Biretta et al. 1996). In addition, the stel- 
lar positions (section 3.3) were determined from 
dithered images taken through the F218W and 
F547M filters. The pixel size was 0.0455 arcsec 



The references in Richards et al. (1999) to a binary sep- 
aration of 25 au are made in error - the binary orbital 
semi-major axis is 25 au, giving the binary separation as 
50 au. 
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Table 1: Observation log. 



Filter 


Exposure 


A 


AA 


Peak A 


Scientific features 




times (s) 


(A) 


(A) 


(A) 


and wavelengths (A) 


F218W 1, 


40 


2136 


355.9 


2091 


Interstellar absorption 


F218W 


100 


2136 


355.9 


2091 


Interstellar absorption 


F437N 


800 


4369 


25.2 


4368 


[O III] A4363 


F469N 


100 


4695 


24.9 


4699 


He II A4686 


F487N 


200 


4865 


25.8 


4863 


H/3 A4861 


F502N 


100 


5012 


26.8 


5009 


[O III] AA4959, 5007 


F547M+ 


2 


5454 


486.6 


5362 


Stromgren y 


F547M 


20 


5454 


486.6 


5362 


Stromgren y 


F656Nt 


100 


6562 


22.0 


6561 


Ha A6563 



t dithered images 



for the undithered images and 0.02275 arcsec for 
the dithered ones. 

The VLA observations were made on 1999 Septem- 
ber 26 at 8.56 GHz and 23 GHz. Comparison 
with the primary calibrator 1331+305 (3C286) 
gave the flux of secondary calibrator 1935+205 
as 0.385+0.001 Jy at 8.56 GHz and 0.41+0.01 Jy 
at 23 GHz. The complex gain solutions for this 
calibrator were applied to HM Sge. MERLIN ob- 
servations by Richards et al. (1999) show that the 
structures in the inner nebula move by ^4 milli- 
arcsec per year and the flux density remained 
roughly constant at 22 GHz between 1994 and 
1996. Thus the VLA and HST images can be 
directly compared. 

3. Results 

3.1. HST WFPC2 images 

The HST WFPC2 images are shown in Figs. 1 
and 2. Some extended emission is apparent out to 
2 arcsec, but the brightest parts are within 0.5 arc- 
sec of the central star, which is readily identified 
from the diffraction spikes in Fig. 1(c). Feature 1 
visible in Fig. 1(a) is almost a ring of emission at 
the south-west edge of the inner nebula. There 
are three compact knots to the east of the central 
emission (features 2, 3 & 4), and a fourth one to 
the south-west (feature 5), which arc most obvi- 
ous in the F502N filter [Fig. 1(a)]. Weak filaments 
are also present to the north (feature 6 & 7), and 
may be associated with northern prominences seen 



at larger scales from the ground by Corradi et al. 
(1999). The inner parts of the nebula shows three 
structures: a ridge to the north (N), a peak at 
the center (C) and a more irregular feature to the 
south (S), which may be the brightest part of a 
loop of emission. In the case of the F656N image 
Fig. 2(a), dominated by Ha, the central star posi- 
tion has been fitted and the point spread function 
(PSF) subtracted, using a PSF modeled using the 
TinyTIM software (Krist 1995). This image shows 
that features C, N and S are distinct in hydrogen, 
and they are also present in the weaker H/3 emis- 
sion in the F487N image Fig. 1(c). The features C, 
N and S are present in the F218W image Fig. 2(d), 
but the latter two features are extremely weak. 

3.2. VLA images 

The VLA images are presented in Figs. 3(a) & 
(b). At 23 GHz there is a northern ridge running 
E-W, a central, almost circular peak, and a more 
irregular southern feature. These features coincide 
with features N, C and S in Fig. 1(a), and we adopt 
the same labels for the radio structure. From com- 
parison of Fig. 3(a) and e.g. Fig. 1(a) we suggest 
that at least one of the binary components is co- 
incident with radio feature C. We discuss this fur- 
ther in section 3.3. The 8.56 GHz image [Fig. 3(b)] 
has significantly lower resolution, but is sensitive 
to more extended emission. Two features are ap- 
parent: a central N-S elongated structure, encom- 
passing features N, C and S, and an extension to 
the S W coincident with feature 1. 
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Fig. 1.— HST images of HM Sge, dereddenned for E(B-V) = 0.35, in the WFPC2 filters (a) F502N; 
greyscale range 9.1xl0~ 15 to 1.3xlCr 12 (main panel) and 1.3xfCT 12 to 2.9X1CT 11 (inset); (b) F437N; 
greyscale range 4.7xlCr 15 to 2.3xfCT 13 (main panel) and 2.3xl0~ 13 to 10" 11 (inset); and (c) F487N; 
greyscale range 3.8xl0 -15 to 9.6xl0~ 13 (main panel) and 9.6xl0~ 13 to 2.5xlO -11 (inset). Units are 
erg s _1 cm -2 A -1 arcsec~ 2 . In each case the inset is at the same scale as used in the radio images in 
Figs. 3(a) and (b). Image (d) is the ratio of images (a) and (b) with greyscale range 2 to 50, and a single 
contour at Rom = 15- Uncertainties in Rqiii are of order 10%. 



4 




Fig. 2.— HST images of HM Sge, dereddenned for E(B-V) = 0.35, in the WFPC2 filters (a) F656N 
star subtracted; greyscale range 1.47xlCr 15 to 2.93xfCT 14 (main panel) and 2.93xl0 -14 to 2.14xl0 _12 
[inset, same scale as Figs. 3(a) and (b)]; (b) F547M; greyscale range 1.58xfCT 16 to 4.74xfCT 13 ; (c) F469N; 
greyscale range 2.09xfCT 16 to 1.47xlO -11 ; and (d) F218W; greyscale range 2xl0~ 13 to 2xf0~ 12 . Units are 
erg s _1 cm~ 2 A -1 arcsec" 2 . In each case the inset is at the same scale as used in the radio images in Fig. 3. 
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Fig. 3.— Radio images of HM Sge at (a) 23 GHz with the VLA; contours are -3, 3, 6, 12, 18, 24 and 36 x 
262 fiJy beam -1 , greyscale range is 0.262 to 12.05 mJy beam -1 . In each case the beam size is illustrated in 
the bottom left hand corner; and (b) 8.56 GHz with the VLA; contours are -3, 3, 6, 12, 24, 48, 96 and 192 x 
80 /iJy beam -1 , greyscale range is 0.08 to 18.96 mJy beam -1 , (c) Extinction map (greyscale) derived from 
Figs. 1(c) and 3(a), with VLA image at 23 GHz overlaid as contours. Propagated uncertainty in E(B-V) 
is ^0.07. (d) Electron temperature versus density; solid lines are the loci for values of Rom as indicated. 
Dashed lines are for v = 8.56 GHz, 7], = 140 K [lowest contour in Fig. 3(b)] and different values of I (i) 
0.00063 pc, (ii) 0.005 pc, (iii) 0.025 pc. Dotted lines are for v = 8.56 GHz, T b = 2260 K [fifth contour in 
Fig. 3(b)] and different values of I (iv) 0.00063 pc, (v) 0.005 pc, (vi) 0.025 pc. The dot-dashed line (vii) is 
for v = 23 GHz, T(, = 6760 K [the peak of feature C in Fig. 3(a)]. See text for discussion. 
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3.3. Stellar positions 

A simple estimate of the relative contributions 
of the HC and CC in HM Sge at a given wavelength 
can be made assuming blackbodies with = 
200 000 K (Miirset, Wolff & Jordan 1997), T^p 
= 3 000 K, and taking approximate radii of 0.01 
and 100 R Q respectively. This shows that the HC 
dominates the CC by a factor ~ 80 at 220 nm 
[F218W, Fig 2(d)] and the CC dominates the HC 
by a factor - 18000 at 550 nm [F547M, Fig 2(b)]. 
Note that UV spectra show negligible line contri- 
bution to the F218W filter (see e.g. Kenyon 1986). 
A number of lines contribute to the F547M filter 
(see e.g. Schmid et al. 2000, their Fig. 2). The 
CC is clearly badly modeled by this estimate, as 
it does not account for the absorption bands char- 
acteristic of Mira-types. In addition, a simple es- 
timate of the luminosity of the HC Miirset et al. 
(1991) and the contribution at A ~2200A shows 
that the central pixel will be entirely dominated 
by stellar emission, rather than nebular emission. 
These estimates show that if the two stars were 
displaced by more than a few pixels on the sky, 
the peak of the emission should clearly move be- 
tween the F218W image and the F547M image. 
That this is not the case means that both stars 
must be within feature C. 

As the PSF of the star in each image extends to 
a large angular distance from the peak, we have 
attempted to fit a Gaussian component to the 
peak positions in the F218W and F547M images. 
The best resolution was available in our dithered 
images, where the pixel size is 0.02275 arcsec. 
The position of the peak relative to 19 41 57 
+16 44 39 in the F218W image [Fig. 2(d)] 
is a = 0.08067±0.00004, 5 = 0.6195±0.0005 
while in the F547M image [Fig. 2(b)] it is a = 
0.08258±0.000005 S = 0.5920±0. 00002. The er- 
rors quoted are the formal uncertainties of the fits. 
When dithering, the measured shifts (from the po- 
sitions of bright stars in each frame) can be up to 
0.2 undithcrcd pixels or ~ 9 milli-arcsec (mas) 
different from the requested value of 5.5 pixels. 
This is probably the best estimate we have of the 
true uncertainty of the measured peak shift from 
F218W to F547M. Thus, the peaks in the two fil- 
ters are displaced from each other by 40±9 mas 
at a position angle of 130°±10° measured north 
through east. The position angle agrees well with 
that suggested by Schmid et al. (2000). We sug- 



gest that the peak position in the F218W image 
is that of the HC while the peak position in the 
F547M image is that of the CC. Thus, the line 
connecting the two stars on the sky, the "binary 
axis", is at a position angle of 130°±10°. The 
schematic in Fig. 4 illustrates the relationship be- 
tween the binary components and the features 
apparent in the inner nebula. 

To test further these measurements, we also 
looked at other observations of symbiotic stars 
made under our GO programme. In the case 
of CH Cyg (paper in preparation), the star was 
known to be in eclipse during the observation, so 
that only one star could possibly be seen. The 
shift in that case was ~ 10 mas, comparable to 
the estimate from the dithering uncertainties, and 
hence consistent with the star being in eclipse. 

The peak positions in the other filters are not 
necessarily coincident with the stars as they are 
dominated by nebular lines. These positions are 
illustrated in Fig. 4. By definition these filters are 
sensitive to more extended emission. This means 
that fitting a gaussian is invalid, as the nebulosity 
may be asymmetrical but strongly emitting near 
the stars. Thus a lack of coincidence between the 
peaks of the narrow filters and the peaks of the 
F218W and F547M filters is to be expected. The 
Ha and H/3 images have peaks very close to that of 
the F547M image. The He II image, which traces 
the highest ionisation regions, has a peak between 
the two stars and slightly south of the binary axis, 
consistent with the model of Nussbaumer & Vo- 
gel (1990). Finally, the F437N and F502N images, 
dominated by [O III] lines which are sensitive to 
temperatures and densities, have peaks which dif- 
fer from those of all the other images and from 
one another. The results for these latter three 
filters are most readily explained in the context 
of the suggested wind interaction which might be 
expected to have the strongest effects between the 
two stars. 

3.4. Nebular diagnostics 

A number of diagnostics of the physical condi- 
tions in the nebula can be derived from the images 
presented here. The ratio Rom of the F502N im- 
age (including both [O III] AA4959 & 5007) to the 
F437N image (including [O III] A4363) depends 
on both electron temperature T e and density n e 
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Fig. 4. — Schematic illustrating the relationship between the binary components and the features of the 
inner nebula (not to scale.) North is up. The position angle of the binary axis is #=130°±10°. The inner 
nebula features (N, C, S and 1) are marked. The hot component (HC) and the cool component (CC) are 
marked at the positions of the peak in the F218W and F547M image respectively. Other peak positions are 
(i) F437N, (ii) F469N (iii) F487N & F656N and (iv) F502N, marked as crosses. 
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according to the equation 

1X0111 ~ JA4363 „ , (U 

_ 7.73exp[Z.2QxW i T~ 1 } \ L ) 
~ l+4.5xlO- 4 (n e T e - 5 ) 

(see Osterbrock 1989). Examination of ground- 
based optical spectra indicate that these two 
WFPC2 filters are dominated by the three lines 
which contribute to the ratio, and the nebular con- 
tinuum is very weak (see e.g. Schmid et al. 2000, 
their Fig. 2). Thus, Figs. 1(a) and (b) can be used 
to derive the ratio map Fig. 1(d) which traces T e 
and n e , and in principle depends on both of these 
quantities. This diagnostic has been applied ef- 
fectively to HM Sge by Schmid & Schild (1990), 
using spatially-unresolved spectra for the entire 
nebula. 

The WFPC2 Charge Transfer Efficiency (CTE) 
effect on the WFPC2 chips (Whitmore & Heyer 
1998; Whitmore 1998) has a bearing on this anal- 
ysis. For aperture photometry, this causes stars 
at the top of the chip (row number 800) to appear 
systematically fainter than those at the bottom 
(row number 0). The effect is also stronger for 
brighter stars. The consequences of this effect for 
extended emission is not well understood. In our 
case the observed emission falls on the same re- 
gion of the chip for all observations, and brightness 
variations are not severe across the nebula. This 
means that we can draw firm qualitative conclu- 
sions from Fig. 1(d). 

Radio brightness, expressed as a brightness 
temperature T&, is given by 

r 6 = T e (l-e- T ") K, (2) 

where the optical depth at frequency v is 

r v ~ 8.24 x 10- 2 T- 135 u- 21 n 2 e l (3) 

assuming n e is constant, where I is the path length 
through the nebula (again see Osterbrock 1989). 
Thus, the radio brightness temperature also de- 
pends on both T e and n e . Relationships (1), (2) 
and (3) provide loci in a T e -n e plot, which inter- 
sect at the conditions appropriate to the various 
features. Comparing Figs. 1(d) and 3(b), we can 
determine the conditions as a function of position. 

Loci for typical values of Rom, Ti> and I are il- 
lustrated in Fig. 3(d). Also shown here is the locus 
for the Schmid & Schild (1990) value of Rqiii = 



5.9. It is clear that the main source of uncertainty 
in this analysis is the value adopted for I. A rea- 
sonable estimate would be to place it at some- 
where between the largest and the smallest angu- 
lar size seen on the sky. When estimated in this 
fashion, it also depends linearly on the distance, 
which is not well known. However, we should still 
be able to draw qualitative conclusions about the 
relative temperatures and densities of the different 
regions despite this distance uncertainty. Most im- 
portantly, it is apparent from Fig. 3(d) that along 
contours of constant Fig. 1(d) primarily traces 
variations in T e . At the same time, variations in Tf, 
trace variations in n e . We must be careful not to 
use this analysis to draw conclusions about those 
regions of the nebula where no radio emission is 
detected. An interferometer such as the VLA will 
resolve-out smooth extended structure, meaning 
that the true extended radio brightness will gener- 
ally be greater than that present in the radio map. 
We note that the value of Rom from Schmid & 
Schild (1990) is effectively an average, dominated 
by the brightest nebular feature. We have plot- 
ted a locus in Fig 3(d) for the peak emission in 
the 23 GHz image (the peak of feature C) as the 
most relevant one for comparison with the Schmid 
& Schild (1990) locus. Given the variable nature 
of the inner nebula (Richards et al. 1999), the val- 
ues which might be derived for that central region 
are similar to those derived by Schmid & Schild 
(1990), which were effectively a weighted average 
for the entire nebula. 

The F469N image [Fig. 2(c)] is dominated by 
He II A4686 emission. This is a discriminator of 
strongly ionized regions (see Osterbrock 1989). It 
can be seen from the image that features C, N and 
S are visible in He II. This suggests that these fea- 
tures are strongly ionized, and is consistent with 
these features being well defined in the F437N and 
F502N images [Figs. 1(a) & (b)], as [O III] shows 
roughly the same ionization structure. In addi- 
tion, the He II emission does not extend far be- 
yond a radius of ~1 arcsec, suggesting that the 
intermediate scale nebulosity is more weakly ion- 
ized than the smallest scale features. This is more 
difficult to reconcile with the extended structure 
seen particularly in Fig. 1(a). However, as noted 
above the [O III] lines seen in these filters are sen- 
sitive to temperature and density, so we may well 
be seeing the effects of density inhomogeneities in 
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the nebula. Such inhomogeneities are in turn con- 
sistent with the very clumpy nature of much of the 
emission at sub-arcsec scales. 

3.5. Extinction mapping 

A third diagnostic comes from the fact that the 
dereddenned H/3 flux, F(H/3), can be derived from 
the radio flux S v via the equation 

S*„ = 2.51 x 10 8 T , c °- 53 ^"°- 1 J F(H/3) Jy, (4) 

where F(Hj3) is in erg cm _2 s _1 . This assumes the 
radio flux is due to optically thin thermal emission. 
Thus, with a measured H/3 flux we can derive the 
reddening E(B-V). In the past this has been used 
to derive the interstellar reddening using the total 
radio and H/3 fluxes. However, in dusty environ- 
ments such as that seen in HM Sge, there may be 
a contribution to the reddening from circumstellar 
matter. Such variations in circumstellar redden- 
ing have been suggested to explain the different 
estimates of E(B-V) and Ay towards HM Sge in 
the past (e.g. Munari & Whitelock 1989). In our 
case, the images of both H/3 emission [Fig. 1(c)] 
and radio emission [e.g. Fig. 3(a)] allow us to de- 
rive an extinction map, as is shown in Fig. 3(c) 
for T e = 10 000 K. We note that Ivison, Hughes & 
Bode (1992) present a spatially-unresolved radio- 
infrared spectrum for HM Sge, indicating that the 
turn-over from optically thick to thin radio emis- 
sion occurs at 8.5 GHz. However, it is clear from 
their Fig. 1 that the turn-over is gradual, con- 
sistent with partially optically thin emission up to 
~30 GHz. In addition, the inner nebula of HM Sge 
is variable in both brightness and structure (e.g. 
Richards et al. 1999), making the spectrum of Ivi- 
son, Hughes & Bode (1992) out of date. This 
means that it remains difficult to assess the op- 
tical depth conditions within the inner nebula as 
a function of position. 

The consequences of the CTE effect for Fig. 3(c) 
are to provide an overestimate of the extinction, 
with the largest over estimates being at the bright- 
est parts of Fig. 1(c). This extinction map is 
also affected by our assumptions about the ra- 
dio emission. Equation 4 relics on the assumption 
that the radio emission is optically thin thermal. 
While the brightness temperatures are consistent 
with thermal emission, the optical depth is not 
well constrained. However, where emission is op- 
tically thick, we would underestimate the extinc- 



tion. As the brightest radio emission most likely 
comes from the regions of greatest optical depth, 
this effect works in the opposite sense to the CTE 
effect. Finally, E(B-V) depends on electron tem- 
perature. If T e = 40 000 K, instead of the 10 000 K 
assumed above, E(B-V) would decrease by <~ 0.2. 
Fig. 3(d) demonstrates that T e can vary consid- 
erably over the nebula. However, for the larger 
values of Tb in the inner regions relationships (2) 
and (3) lead to loci at higher electron densities, 
while the curves derived from equation (1) con- 
verge towards lower T e , and extend over a nar- 
rower range of temperatures. Thus in the inner 
regions shown in Fig. 3(c) the variations in T e are 
relatively small. As the hotter regions are closest 
to the white dwarf, it seems likely that the inter- 
mediate values of E(B-V) = 0.7 may be somewhat 
reduced by the expected higher T e in that region. 
Bearing in mind these considerations, we are confi- 
dent that we can draw qualitative conclusions from 
Fig. 3(c). 

4. Discussion 

4.1. Dust distribution and interstellar ex- 
tinction 

The extinction map in Fig. 3(c) can be taken to 
trace the circumstellar variation in the dust dis- 
tribution, assuming the interstellar extinction is 
reasonably uniform over the relatively small an- 
gular size of the inner nebula. In this case, the 
best estimate of the true interstellar extinction 
comes from the minimum values in the extinction 
map. Three distinct regions are apparent: (i) to 
the south E(B-V) ~ 1, coincident with feature S; 

(ii) to the north around feature N, E(B-V) ~ 0.6, 
consistent with previous reddening estimates; and 

(iii) a low-ratio band running from east to west 
across the position of feature C. Note that the 
peak itself does not give a reliable E(B-V) esti- 
mate, as the emission in the optical is not entirely 
nebular. This band is also weakly contaminated 
by the diffraction spikes, which have a position 
angle of ~ 80°. The eastern part of this low- 
ratio band has E(B-V) = 0.35, the lowest value 
on the map, and presumably the best estimate of 
the true interstellar extinction. As high E(B-V) 
correlates with a greater quantity of dust, we can 
see that the dust is concentrated at feature S. In 
addition, there is a deficit of dust closer to fea- 
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ture C in Fig. 1(a), especially immediately to the 
east. This clumping of the dust at small scales was 
suggested by Thronson & Harvey (1981), based on 
discrepancies in reddening estimates. Our maxi- 
mum value of E(B-V) ~ 1 gives Ay ~ 3.1 for in- 
terstellar dust (Evans 1994), which is much lower 
than the value of Ay ~ 12 of Thronson & Harvey 
(1981) but similar to other estimates both before 
and after that time. Much more recently, spectro- 
polarimetry by Schmid ct al. (2000) demonstrated 
that the Mira was clearly visible, in agreement 
with the significant reduction in Ay since 1980 we 
see here. Results from the Infrared Space Observa- 
tory (Schild et al., submitted) are consistent with 
a two-component dust model, with optical depths 
comparable to those derivable from our E(B-V) 
values. Such a model is also consistent with the ex- 
tinction map in Fig 3(c), particularly as the Schild 
ct al. modeling suggests that at least one of the 
components is clumpy. 

We suggest that the highest concentration of 
dust, at feature S, is associated with the CC wind. 
Note that Fig. 3(c) is blanked beyond the la level 
of the VLA image, as no useful information can 
be derived beyond this limit. Hence, we have 
no knowledge about the dust distribution beyond 
that region. 

4.2. The distance and the binary parame- 
ters 

The minimum in the extinction map Fig. 3(c) 
places the interstellar extinction at no more than 
E(B-V) = 0.35. This is an upper limit because 
there may be a circumstellar contribution to the 
minimum. This gives an interstellar total extinc- 
tion A v <1.1 magnitudes. Neckel & Klare (1980) 
provide maps of the galactic extinction versus dis- 
tance. HM Sge lies at the edge of their field 265 
(see Figs. 5 and 6m of Neckel & Klare 1980). For 
Ay = 1.1, the minimum distance to HM Sge is 
approximately 700 pc, although measurements of 
similar Ay are shown in field 265 for distances as 
low as 400 pc and as high as 2 kpc. Comparison 
with the extinction maps of Lucke (1978) suggest a 
distance of ~800 pc for this value of E(B-V). How- 
ever, such large-scale extinction maps are based 
on a small number of lines of sight in each field, 
so that the distance-dependent extinction values 
are not readily applicable to individual objects for 
the purposes of distance estimates. 



Our measured binary angular separation (sec- 
tion 3.3) suggests a projected binary separation of 
(40±9)D au at a distance of D kpc. This can be 
reconciled with the results of Richards et al. (1999) 
if the distance to HM Sge is ~ 1250±280 pc, which 
is consistent with our estimates above. 

4.3. Nebular temperatures 

The map of Rom [Fig. 1(d)] shows two dis- 
tinct regions in the extended nebula visible with 
the HST. Most of the nebula has 4 < Rom < 
15, while a wedge region to the south-west has 
Rom > 15 [coincident with feature 1 in Fig. 1(a)]. 
At the center, where the stars contribute signifi- 
cantly, Ron i is as low as 2, but clearly this is not 
purely nebular emission, and so cannot be relied 
upon for nebular diagnostics. From Fig. 3(d) it 
is clear that higher Roin regions are consistent 
with lower electron temperatures along a given Tb 
contour. This wedge extends from roughly the lo- 
cation of the highest extinction in Fig. 3(c) to the 
edge of the nebula, and is suggestive of the "bow 
wave" structures seen in interacting-wind simula- 
tions. More importantly, there is a cooler region 
in the nebula associated with the southern dust 
feature, further supporting our suggestion that we 
are seeing the CC wind directly. Alternatively, 
this region of the nebula could be shielded from 
the HC radiation field by feature S. 

5. Conclusion 

Using contemporaneous HST WFPC2 and 
VLA observations of the symbiotic nova HM Sge, 
we have investigated the nebular conditions and 
the relationship of the star positions to the neb- 
ular structure. We identify a number of discreet 
features at spatial scales smaller than ~ 0.1 arcsec 
embedded in the extended nebula, including at 
least four discreet knots, and two linear features 
to the north-east. A ring-like feature (labeled 1) 
is evident on the south-west edge of the inner 
nebula. 

The radio and optical emission is well correlated 
in the inner 1 arcsec, with a northern ridge (N), 
a central peak (C), and a more irregular southern 
feature (S) clearly evident at both 23 GHz and 
^5000A. Temperature diagnostics show a "wedge" 
to the south-west, encompassing feature 1, con- 
sistent with lower temperatures than the rest of 
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the nebula. The knots 1-3 and the western tip of 
ridge N also appear to have cooler, denser condi- 
tions than the rest of the nebula, but these are less 
well defined in the diagnostic map. 

Both stars must be encompassed by fea- 
ture C. The peak shifts by 40±9 mas between the 
F218W image (~ 220 nm) and the F547M image 
(~ 550 nm), with a position angle of 130°±10°. 
We place the HC at the continuum UV peak and 
the CC at the red continuum peak. The posi- 
tion angle of the line between the two stars on 
the sky, the "binary axis", is 130°±10°. This 
is consistent with a suggested binary axis from 
spectro-polarimetry (Schmid et al. 2000). It is 
also consistent with a suggested binary separation 
of 50 au (Richards ct al. 1999) if the distance is 
- 1250±280 pc. 

The circumstellar extinction map suggests a 
patchy dust distribution, with the greatest extinc- 
tion to the south-west, coincident with feature S, 
and the "wedge" around feature 1. The mini- 
mum in the extinction map indicates the true in- 
terstellar extinction to be no more than E(B-V) = 
0.35. This is consistent with a minimum distance 
of ~700 pc, but this would be reduced if there 
is a circumstellar contribution to the minimum. 
We suggest feature S and the south-west wedge 
(feature 1) are associated with the CC wind. Al- 
ternatively, the south-west wedge is shielded from 
the HC radiation field by feature S. 
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